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ABSTRACT 

We used ultra-deep UV observations obtained with the Hubble Space Tele¬ 
scope to search for optical companions to binary millisecond pulsars (MSPs) in 
the globular cluster 47 Tucanae. We identihed four new counterparts (to MSPs 
47TucQ, 47TucS, 47TucT and 47TucY) and conhrmed those already known (to 
MSPs 47TucU and 47TucW). In the color magnitude diagram, the detected com¬ 
panions are located in a region between the main sequence and the CO white 
dwarf cooling sequences, consistent with the cooling tracks of He white dwarfs 
of mass between 0.15 Mq and 0.20 Mq. For each identihed companion, mass, 
cooling age, temperature and pulsar mass (as a function of the inclination an¬ 
gle) have been derived and discussed. For 47TucU we also found that the past 
accretion history likely proceeded in a sub-Eddington rate. The companion to 
the redback 47TucW is conhrmed to be a non degenerate star, with properties 
particularly similar to those observed for black widow systems. Two stars have 
been identihed within the 2a astrometric uncertainty from the radio positions 
of 47TucH and 47TucI, but the available data prevent us from hrmly assessing 
whether they are the true companions of these two MSPs. 

Subject headings: Pulsars: Individual: J0024—7204H, J0024—72041, J0024—7204Q, 
J0024-7204S, J0024-7204T, J0024-7203U, J0024-7204Y, Globular clusters: 
Individual: 47 Tucanae (NGC 104), Techniques: photometric 


^Based on observations collected with the NASA/ESA HST (Prop. 12950), obtained at the Space Tele¬ 
scope Science Institute, which is operated by AURA, Inc., under NASA contract NAS5-26555. 
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INTRODUCTION 


Millisecond pulsars (MSPs) are rapidly spinning neutron stars (NSs) formed in a bi¬ 
nary system where a slowly rotating NS is spun up through mass accretion from an evolv¬ 
ing companion star. The recycling process is usually observed in the l ow mass X-ray bi¬ 


nary systems, which are commonly considered as the M SP progenitors flAlpar et ahlll982 


Bhattacharva fc van den Heuvellll99ll:lPapitto et al.ll2013[l . Indeed, when the mass accretion 
rate decreases and the NS is sufficiently recycled, the rotation-powered emission is switched 
on and the system is observable as a MSP in the radio band. These processes usually lead to 
a deep transformation of the companion star which ca n transit through a highly perturbed 
evolutionary phase (possibly like MSP-A in NGC 6397: iFerraro et al.ll2001al l. b efore reaching 
the fin al stage of a (possibly He) white dwarf (WD; e.g. MSP-A in NGC 6752: iFerraro et ah 
2003all . Each step of this evolution corresponds to objects characterized by different prop¬ 
erties, that, in the optical bands, are imprinted in the observable features of the companion 
star. 


According to the can onical scenario, the majo rity of binary MSPs have low-mass He 
WD companions (see, e.g.. Ivan Kerkwiik et al.ll200511 . However, recent pulsar (PSR) searches 
have considerably increased the number of non-canonical systems, especially the so-called 
“black widows” and “redbacks”: ultra-compact binary systems (Pqrb 1 day) where the 
presence of radio eclipses suggests the presence of ionized material ablated from a bloated 


companion star beca use of the energy injected by the PSR flRuderman et al.lIlQSQuRav et ah 


2 OI 2 I: lRobertsll2013h . Redback companion stars usually have masses of 0.1 — 0.5 Mq, while 


black widow companions are much less massive (M <0.1 Mq). Such a small value could 
be due to vaporization from the strong MSP radiation and relativistic wind. The physical 


mecha nisms bringing to one or the other system are still debated. Simulations bv iGhen et ah 


(I 2 OI 3 I) show that redbacks and black widows are the outcome of different evolutionary paths, 
where the PSR irr adiation efficiency is th e discriminant factor. At odds with these results, 
the simulations bv iBenvenuto et ahl (120141 ) show that the evolution of redbacks is bifurcated, 
with some of them evolving into black widows, and the others producing canonical He WD 
systems. Possibly, the progressive evaporation of the black widow companions could lead to 
the total disruption of the star and then to the formation of isolated MSPs. Interestingly, 
in recent years several connections between low mass X-ray binaries and redbacks have been 


Archibald et ah 

2009; 

Panitto et ah 

2013; 

Bassa et ah 

2014) 


Although the Galaxy is ~ 100 times more massive than the entire Galactic glob¬ 
ular cluster (GG) system, about 40% of the known MSP population is found in GGs. 
Such an over-abundance is indicative of a strongly enhanced dynamical activity in these 
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dense stellar systems, which promotes the formation of a conspicuons nnmber of ex otic ob¬ 
jects, such as blue straggler stars. X-ray binaries , cataclysmic variables and MSPs flBailvn 


1992 


2003 


Cool et af 


1995 


Ransom et ah 


Ferraro et al.l Il995l. l2001bl l2015ai ICrindlav et ahl l2002l: IPoolev et al 


20051) ■ which can be used to probe the complex i nterp l ay between dy¬ 


namics and stellar e volution (e.g. 


Possenti et ah 2003: Ferraro et ah 


Goodman fc Hut 


2003b 


2009 


11989 


20121 


Hut et ah 


2015b 


1992 


Phinnevl 1992 


Verbunt fc Freirel 12014) . In 


this respect, the study of optical companions to binary MSPs in GCs is of outmost im¬ 
portance, since it opens the possibility to get insights on the impact of dynamical inter¬ 
actions (which are particularly frequent in dense environments) on MSP and stellar evo¬ 
lution, e.g., favoring binary formation (through tidal captures), binary shrinking (through 
fly-by) and consequent mass transfer activity, as well as exchange interactions able to sub¬ 


stitute the original companion t 

lat recvciec 

the pulsar, with a new, rr 

star (see e.g. 

Rasio et al 

2000l: 

King et ah 

2003 

Ferraro et al. 

2003c 

Mucciarelli et ah 

20131 E 

lenacQuista & Downingl 

2 OI 3 I). Moreover, in 


panions, it is possible to estimate the masses and cooling ages of the s ystems by the direct 
comparison of their p roperties with stellar evolutionary models (see e.g. iFerraro et al.ll2003al: 
Pallanca et al.ll2013al). while accurate mass measurements require spectroscop ical techniques 


( e.g. Ivan Kerkwijk et al.lll996l:lBassa et al.ll2nn6l: lAntoniadis et al.ll2n 121 . 120131 ). In the case of 
MSPs in GCs, these constraints benefit from the known GC distance and optical extinction, 
thus reducing the uncertainties on the estimated quantities with respect to the case of MSPs 
in the Galactic field. The derived companion masses can be combined with radi o timing pa¬ 


rame ters to estimate the PSR masses, thus allowing general relativity test (e.g. iFreire et al. 


20121) and fun damental physics studies, a s the determination of the equation of state of ultra 
dense matter (ILattimer fc Prakashll200ll) . On the other hand, the derived cooling ages are a 
more appropriate measurer nent of the system age with respec t to the characteristic PSR ages 
derived from radio timing flTauris et al.ll2012al: lTaurisll2012bl ). The correct determination of 
MSP ages is an important to ol to study the spin evolu tion and to constrain the physics of 


the recycling phases (see e.g. Ivan Kerkwiik et al.ll2005l. and references therein). 


Despite their importance, the identihcation of MSP optical companions is challenging 
in crowded stellar systems like GCs. Only te n companions have been discovered so far 


Sigurdsson et al. 

2003; 

Bassa et al. 

2003, 

2004' 

, as expected 

from 

the canonical 1 

orma- 

tion scenario, five are redbacks companions (see 

Ferraro et al. 

2001a: 

Edmonds et al. 

2002; 


Cocozza et al. 


2008 


(iPallanca et al. 


Pallanca et al. 


2014al: 


2010 


Cadelano et al. 


20151) . 


2013bl) and two are black widow companions 


The GC 47 Tucanae, located at a dis tance of about 4.5 kpc from the Sun, hosts the 
largest population of MSPs after Terzan 5 (IFreire et al.ll2003l: iRansom et al.l 120051) . Indeed, 
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23 radio MSPs have been discovered so far, 14 of which are located in binary system^. Here 
we report the identihcation and the properties of fonr new MSP companions in 47 Tncanae 
and we present the follow-np stndy of two previonsly known companions. In Table [T] we 
report the main radio timing properties of the analyzed objects, which are nsefnl in the 
following discnssions. All the identihed companions are likely canonical MSPs, except one 
which is a redback system. 

In Section 2 we present the nsed photometric dataset and the identihcation of the MSP 
optical connterparts. In Section 3 we discnss in detail the properties of each companion. 
Finally, in Section 4, we snmmarize onr resnlts. 


2. OPTICAL PHOTOMETRY OF THE STAR CLUSTER 


2.1. Observations and data analysis 


In the present stndy the identihcation of the MSP companions has been performed 
throngh an nltra-deep, high resolntion, photometric dataset acqnired nnder GO 12950 (P.I; 
Heinke) with the UVIS camera of the Wide Field Camera 3 (WFC3) monnted on the Hnbble 
Space Telescope (HST). The dataset consists of 8 images in the F390W hlter, with exposnre 
times of 567 — 590 s, and 24 images in the LP F300X hlter, with exposure times of 604 — 609 s. 


The standard photometric analysis (see iDalessandro et al.ll2008alJbl ) has been performed 
on the “ht” images, which are corrected for hat held, bias and dark counts. These images 
have been further corrected for “Pixel-Area-Map”§ with standard IRAF procedures. By 


ap L 

using the DAOPHOT II packages (iStetsonl 119871) . we performed an accurate photometric 
analysis of each image. First of all, we modeled a spatially varying Point Spread Function 
(PSF) by using a sample of ~ 200 bright but not saturated stars. The model has been chosen 
on the basis of a test and, in every image, the best ht is provided by a Mohat function 
fjMohatlll969l) . Then we performed a source detection analysis, setting a 3cr detection limit, 
where a is the standard deviation of the measured background. Once a list of stars was 
obtained, we performed a PSF-htting in each image using the ALLSTAR routine. In the 
resulting catalog we included only objects present at least in half the images for each hlter. 
Then, this catalog has been further processed with the ALLFRAME routine. For each star, 
we homogenized the magnitudes estimated in diherent images, and their weighted mean and 


^Please visit http;//www.naic.edu/~pfreire/GCpsr.html, for a complete list of the main radio timing 
properties of MSPs in GCs. 

^For more details see the WFC3 Data Handbook. 
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standard deviation ha ve been finally adopted a s the star mean magnitnde and its related 


photometric error (see iFerraro et ahl Il991 


1992I) . However, in order to perform variability 


studies, for each source we also kept the homogenized magnitude measured in each frame in 
both hlters. Then, instrumental magnitudes have been calibrated to the VEGAMAG system 
by using the zero points quoted in the WFG3 Data Handbook and by performing aperture 
corrections. 


2.2. Astrometry 


Since the WFG3 im ages suffer from geom etric distortions, we corrected the instrumental 
positions (x,y) following iBellini et ahl (1201 11 1. In order to transform the instrumental posi¬ 
tions into th e absolute astr o metri c system (a, 5), we used, hrst of all, the wide held catalog 
presented in iFerraro et ahl (1200411 . Its astrometric solution has been improved b y cross- 
correlatioE0 with the UGAG4 astrometric standard catalog f Zacharias et ah 2013 : ~ 4600 
stars have been found in common between the two datasets). The latter is based on the 
International Gelestial Reference System, thus allowing a more appropriate comparison with 
the MSP positions derived from timing using solar system ephemerides (which are referenced 
to the same system). The newly-astrometrized wide held catalog has then been used as a 
secondary reference frame to astrometrize the WFG3 data set, by means of ~ 22000 stars 
in common. The resulting la astrometric uncertainty is 0.10" and 0.11" in a and 6 , respec¬ 
tively. Thus the hnal total astrometric uncertainty is ~ 0.15". Unfortunately, there are only 
few stars in common between the WFG3 and the UGAG4 catalogs, since the latter does not 
cover the cluster central regions. This prevented a direct cross-correlation between the two 
catalogs and thus we could not take into account the stellar proper motions between the two 
observation epochs, which would have reduced the astrometric uncertainty. 


2.3. Identification of the MSP companions 


First of all, in order to search for the companions to the MSPs in 47 Tucanae, we 
checked the precision of our ast rometric solution re-identify ing the two companion stars 
already known in the cluster (see iFdmonds et ahl l200ll. 120021 1 . To this aim, we performed a 


^We used CataXcorr, a code which is specifically developed to perform accurate astrometric solutions. 
It has been developed by P. Montegriffo at INAF- Osservatorio Astronomico di Bologna. This package is 
available at http://davide2.bo.astro.it/^paolo/Main/CataPack.html, and has been successfully used in a 
large number of papers by our group in the past 10 years. 
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detailed analysis of all the detectable objects within a 5" x 5" wide region centered on the 
nominal position of each MSP. The companion to 47TucU (COM-47TncU; hereafter all the 
companions will be named as COM-47Tnc followed by the letter of the respective MSP) and 
47TncW have been re-identi£ed in stellar sonrces located at 0.06" from the MSP nominal 
positions. Both the identihcations turn out to be largely within our astrometric uncertainty, 
thus conhrming the accuracy of the adopted astrometric solution. The hnding charts of these 
two reference objects are shown in Figure [H 


Following the same procedure, we se arched for the companions to all the other MSPs 
with a known position flFreire et al.ll2003l. Freire et ah 2015, in preparation). Stars located 
within the 2 a uncertainty from the pulsar position have been considered as possibile counter¬ 
parts. Four companions (to 47TucQ, 47TucS, 47TucT and 47TucY) have been identified on 
the basis of their positional coincidence (all of them are located at a distance < 0.06" from 
the nominal radio position) and of their position in the color magnitude diagram (CMD). 
Two faint stars have been detected also within the 2a uncertainty circle from 47TucI and 
47TucH. However, their distances (0.15" and 0.24" respectively) from the pulsar radio po¬ 
sitions are signihcantly larger than in all the other cases, thus casting doubts about these 
objects being the true optical counterparts (see more discussion in Section The hnding 
charts of all these objects are shown in Figure [Hand their main photometric properties are 
reported in Table [2l Their location in the cluster CMD is shown in Figure [2l where only 
the stars with a sharpness parametei@ |sh| < 0.05 are plotted. As can be seen, with the 
exception of the candidate companion to 47TucI, all the newly identified counterparts are 
located in the region where He WDs are expected, although the candidate companion to 
47TucH could be compatible also with the CO WD cooling sequence (see Section l3^ . Since 
the radio timing properties suggest that these systems are the product of the canonical recy¬ 
cling scenario, their location along the He WD cooling sequences guarantees their connection 
with the MSPs. Note in fact that the probability of a chance coincidence with another He 
WD is extremely low (P ^ 0.1%)|^, since these objects can only be the product of the late 
stage of the evolution of exotic objects like, for example, MSPs and cataclysmic variables. 
The candidate companion to 47TucI is instead a main sequence-like object, and its proper¬ 
ties will be briehy discussed in Section 13.21 As concerns the previously known companions. 


®The sharpness parameter is a DAOPHOT II output that quantifies the stellar-like structure of each 
object fitted with the PSF model. See the User Manual for more details. 

®The chance coincidence probability has been evaluated predicting the number He WD expected within 
a radius equal to the 2cr astrometric uncertainty. We derived the He WD density by direct counting of the 
all objects located among the cooling tracks (see Section ITT] and Figure [3]) and dividing this number by the 
size of the WFC3 field of view. Please note that even including all the stars of the catalog with sharpness 
|sh| > 0.05, the chance probability remains < 0.5%. 
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COM-47TucU is also located along the He WD sequence, while the redback COM-47TucW 
is located in an anomalous region between the main sequence and the WD cooling sequence 
(see Section [231) • 


With the exception of the COM-47TucW, no signihcant variability related to the orbital 
period has been detected. For 47TucU, 47TucY and 47TucW (see Section I3.4p the observa¬ 
tions sample a signihcant fraction of the orbital period. Instead, for the other systems (with 
orbital periods longer than 1 day) the coverage is too poor to allow any appropriate variability 


ate companions 


e.g. 

Stanners et al. 

1999; 

Edmonds et al. 

2002 ; 

Revnolds et al. 

2007; 

Romani & Shaw 

2911: 

Pallanca et al. 

I2ni4a: 

Cadelano et al. 

201 

5), is 


not expected and usually not observed for degenerate objects, since the hux enhancement 
due to re-heating of the companion star by the PSR emitted energy is negligible. 


DISCUSSION 


3.1. The physical properties of the He WD companions 


In order to constrain the main properties of the He WD companions, we have com¬ 
pared the po sition of each c andid ate in the CMD with a set of He WD cooling tracks 
computed by lAlthaus et ah! (120131 ). These models span a mass range from 0.15 Mq to 
0.43 Mq, spaced at about 0.005 Mq for masses between 0.15 Mq to 0.19 Mq and up to 
0.07 Mq for larger masses. We transformed the theoretical luminosities and temperatures 
into the absolute F300X and F390 W magnitudes, by applying the bolometric cop ections 
kindly provided by P. Bergeron fsee lHolberg fc Bergeronll2006ldBergeron et al.ll201ll) . Then, 


the model absolute magnitudes have been tran sformed into 


distance modulus (m — M )n = 13.32 ± 0.10 fiFerraro et al 


E(B — V) = 0.04 ± 0.02 ([Ferraro et ahl Il999l: IZoccali et al 


le ap 


1999 


I 


2001 


arent ones by using the 
and the color excess 


Salaris et al. 


extinction coeffi cients HfsooxMv = 1.77309, Ap^gow/= 1.42879 (ICardelli et al. 


2007) and 


1989 


0’Donnelllll9941) . Figure [3] shows the zoomed portion of the CMD in the WD region with 
a sample of cooling tracks for different masses overplotted. As can be seen, the range in 
mass of the models is large enough to properly sample the portion of the CMD where all the 
companions are located. Therefore we used this set of models to derive the combinations 


Many literatu re works reported on different values of 47 Tucanae distance modulus (see, e.g. 
Woodley et al. and references therein). However, all these possibile values have only a minimal influ¬ 

ence on our derived companion properties (e.g. the derived companion masses would vary of less than ~ 7% 
for all the companion stars) 
























































































of parameters (mass, cooling age and temperatnre) that simnltaneonsly satisfy the observed 
photometric magnitndes in both the hlters, also taking into acconnt the nncertainties on the 
companion magnitndes, distance modnlns and reddening. The best valnes have been evaln- 
ated with a simple statistic. In doing this, linear interpolations (for different masses bnt 
eqnal ages) among the tracks have been performed in order to have a tighter mass sampling. 
We assnmed that each companion is located at the distance of 47 Tncana^ and it is affected 
by the same extinctioij^. Fignre H] shows, for each system, the combination of cooling age 
(left panel), temperatnre (central panel) and PSR mass (right panel) appropriate for the 
derived valne of the companion mass. In particnlar, in each plot the right panel shows the 
resnlts obtained for different valnes of the inclination angle and interesting constrains on 
each system can be drawn. For instance, by setting the inclination angle to 90°, the maxi- 
mnm PSR mass allowed from the inferred companion mass can be evalnated. Conversely, by 
assn ming the minimnm PSR mass eqnal to 1.17 Mq (the lowest mass ever measnred for a 
NS; I Janssen et ahl 120081) . a conservative lower limit to the inclination angle can be derived. 
All these resnlts are also snmmarized in Table [3] where the qnoted nncertainties are the range 
of possibile valnes allowed by the comparison with the theoretical trackJ^. Note that we are 
not analyzing here the cases of 47TncH and 47TncI, which we will discnss in Section [221 


As can be seen, all the companions have masses between ~ 0.15 Mq and ~ 0.2 Mq. 
The derived ranges of ages are in agreement with the lower limits to the PSR characteristic 
ages reported in Table [H The only exception is COM-47TncU, which is discnssed below. In 
principle the mass of COM-47TucQ conld be smaller than onr best-£t valne (0.15Mq), since 
not theoretical tracks for masses below this valne are available. However, already a 0.15 Mq 
companion wonld imply an extremely low valne of the PSR mass (< 1 Mq). This pnzzling 
resnlt conld be partially explained with the difhcnlty of accnrately determine the color of 
the optical connterpart, becanse of the presence of a very close bright object (see Fignre [T]). 

Onr resnlts rnle ont a massive NS in the case of 47TncQ and 47TncT, while the possibility 
of a ~ 2 Mq ns remains opened in the cases of 47TncS, 47TncU and 47TncY. However, 
FignrelHshows that the PSR mass can be signihcantly rednced by assnming an intermediate- 
low inclination angle of the orbital plane. In any case these systems are worthy of fntnre, 
especially spectroscopical, investigations. 


®Even though iFreire et al.l (j200lll measured distance offsets between the cluster MSPs, such differences 
are very small and can be neglected for our goals. 

® The effects of differential reddening are negligible for our goals (see Milone et al. 2012al f3l. 


^°The reader should be aware that the WD parameters should not be assumed at face value as perfectly 
correct but as estimations, since they are model dependent and could also suffer from some hardly quantifiable 
uncertainty linked to the bolometric corrections. 
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We also compared our results with the theoretical predictio ns on the behavio r of th e 
orbital period as a function of the companion mass discussed in iTauris fc Savoniid fll999h. 


Bassa et ah 


Such a model has been already empirically verihed by Corongiu et ah f 2012h and 
(120061) . As shown in Figure El where we also added the two He WD co mpanions identihed 


in NGC 6752 and M4 (see iFerraro et ahl l2003al: ISigurdsson et ahl 120031) . our results are in 


reasonable agreement with the model. The analytical prediction seems to slightly overes¬ 
timate the companion mass or to underestimate the system orbital period. However, this 
model is valid for binary systems with 0.18 Mq < Mwd < 0.45 Mq, thus only marginally 
representative of our sample, where mo st of the companion s appear to be less massive than 


0.18 Mq. More updated models (from llstrate et al 


2014J: gray points in Figure [6]) are in 


better agreement, although they are the results of simulat ions of donor stars w ith metallicity 
Z = 0.02, larger than that of 47 Tucanae (i.e. Z = 0.008, Lapenna et al. 2015 ). 


The brightness of COM-47TucU allowed us to put tighter constraints to the system 
parameters with respect to the other objects. Both it s mass and temperature are in excel¬ 
lent agreement with those reported in iFdmonds et ahl (120011) . while our derived age (~ 0.9 
Gyr) turns out to be 0.3 Gyr larger than their estimate. Such a discrepan cy could be due 
to the different theoretical models used. However, as already noticed by lEdmonds et al. 


( 2001 ). the cooling age is signihcantly lower than the characteristic age of 2.5 Gyrfl This 
discrepancy should not alarm, since the PSR characteristic ages are based on many assump- 
tions and large dev i ation s from the conipanion cooling ages a re commonly observed (see e.g. 
Lorimer &: Kramerl l2012l: iTauris et ahl l2012al: iTaurisI l2012b|) . Using the WD age together 
with the intrinsic spin period derivative (P = 2.7 ± 0.5 x 10“^*^; Freire et ah, in prepara¬ 
tion) and the actual spin period (P ~ 4.343 ms), we evaluated a MSP birth spin period 
(the so-called equilibrium spin period) of Pq ~ 3.576 ms. This value, combined with the 
surface magnetic held (B 3.145 x 10® G) and assuming a NS with a radius of 10 km and a 
canonical mass of 1.4 Mq, can be used to infer the typical accretion rate th at reaccelerated 
the N S during the low mass X-ray binary phase. By using equation (8) of Ivan den Heuvel 
(120091) . we hnd that the system past accretion history likely proceeded at a sub-Eddington 
rate {M/Medd ^ 0.02), as expected frorn the typical evoluti on of close binary systems with 
light donor stars (ITauris fc Savoniielll999l: llstrate et al.ll2014l) . Although the mass accretion 
rate strongly depends on the NS radius, the general result does not change assuming different 
radii or even different NS masses. 


^^This value is based on the estimate of the PSR spin-down rate from the orbital period derivative, which 
is precise enough for this system (Freire et ah, in preparation). 
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3.2. Possible candidate companion stars 


As can be seen from Fignres [2], |3] and [6l the possible companion to 47TucH appears 
to have properties qnite different from those observed for the other companions, hrst of 
all its mnch larger distance from the MSP nominal position (0.24"), which corresponds to 
almost twice onr astrometric nncertainty. Moreover, following the procednre adopted in 
the previous section, we derived for this object a mass of 0.37 ± 0.05 M r:^. This value, 
combined with the binary system total mass of 1.61 Mq flFreire et al.ll2003h . would imply a 
PSR mass of ~ 1.25 Mq, a value slightly lower than expected for a recycled PSR, although 
still acceptable within the uncertainties. Its position in the CMD is compatible also with the 
CO WD cooling sequence, which would increase the probability of a chance coincidence to 
~ 2 — 3%. Furthermore, at odds with the others objects, this candidate counterpart occupies 
an anomalous region in the orbital period companion mass plane shown in Figure |6l Although 
this anomaly could be real (since 47TucH has a large eccentricity, probably due to some kind 
of dynamical interaction), all these pieces of evidence suggest that the observed object is 
probably an isolated WD and the true companion star is still under the detection threshold 
(see Section [33D . 


A possible candidate companion to MSP 47TucI has been also detected (see Figure [T]for 
the finding chart). This is a binary system with a short orbital period (~ 0.23 days) and a 
very small eccentricity. From the PSR mass function, the companion is expected to be a very 
low mass star (Mqom ~ 0.015 Mq). The absence of radio eclipses, probably due to a low 
inclination angle, prevents its characterization as a black widow system. At 0.15" from the 
PSR position, we identihed a star located at the faint-end of the cluster main sequence (see 
Figure [2] and Table [2]). If we assume that the companion is a bloated star seen in a binary 
system with a low inclination angle, such a CMD position could be reasonable. However, 
the lack of any significant variability related to the orbital period prevents us from firmly 
associating this candidate to the MSP. In fact, the orbital period coverage of the F390W 
images is too poor, while the signal to noise ratio of the F300X data allows us to only infer 
that, in case of photometric variability, the maximum variation amplitude must be smaller 
than ~ 0.8 mag. We therefore conclude that it is more likely that the real companion star 
is still under the detection threshold. Indeed, the probability of a chance coincidence with a 
main sequence star is non negligible (~ 45 — 50%). We finally note that another object lies 
within the astrometric uncertainty circle, but its association with the MSP can be excluded, 
since it is a common CO WD, with properties incompatible with the MSP timing ephemeris. 
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3.3. Non detections 


No interesting counterparts have been identified for all the other known binary MSPs. 
These non-detections are likely due to companion stars still under the detection threshold 
(as in the case, e.g., of 47TucE and the black widow 47TucJ), or to the severe crowding 
conditions of the area surrounding the MSP positions (as in the case of the black widows 
47TucO and 47TucR). No search could be performed for 47TucX since its position is outside 
the field of view. 


Considering that the companion to 47TucJ should be a non-degenerate object, its non¬ 
detection in UV passbands cannot be used to get useful information on its properties. In¬ 
stead, the counterpart to 47TucE is expected to be He WD, which remains undetected down 
to our limiting magnitudes (~ 25 in the F300X filter and ~ 25.5 in the F390W filter). Hence, 
taking into account that the cooling age of a ~ 0.17Mm WD at these detection thres holds 
is larger than the cluster age (~ 10 — 11 Gyr; iGratton et al.l 120031: iHansen et al.l 120131) . it is 
unlikely that this star has a mass similar to that estimated for the other companions. It is 
more probable that it is more massive than O. 2 M 0 (corresponding to a faster cooling) and 
its cooling age is larger than 1 Gyr. The same should apply also to the case of 47TucH if 
its true companion is still under our detection limits (as sugge sted above). Interestingly, 
according to the theoretical relation of lTauris fc Savoniiel (1199911 and the orbital periods of 
47TucE and 47TucH (~ 2.3 and ~ 2.4 days, respectively), the companions to both these 
MPSs are indeed expected to have masses > O. 2 M 0 . 


3.4. The companion to the redback 47TucW 


47TucW is the only redback identified, so far, in 47 T ucanae. It is a bina ry MSP with 
a spin period of 2.35 ms, an orbital period of ~ 3.2 hr (IGamilo et al.l 1200011 and a com- 
pa nion mass of ^ 0.15 M ^n. The first optical identification of this system was presented 


m 


Edmonds et al.l (120021) . who suggested that the companion is a perturbed and non de¬ 
generate star with a light curve structure indicating a strong heating by the PSR flux. In 
Figure[7]we show, for both the filters, the light curves we obtained by folding our photometric 
measurements with the most updated radio timing ephemeris (Freire et ah, in preparation). 
The zero orbital phase has been se t at the PSR ascending node timJ^. As ca n be seen, 
in agreement with previous works (iFdmonds et al.l 120021: iBoedanov et al.l 1200511 , the light 
curves present a single maximum-minimum structure, likely due to the heating by the PSR 


^^Please note that we are using a different formalism with respect to lEdmonds et all (l2002fl . 
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flux. Unfortunately, the star has been measured above the detection threshold only near 
its maximum luminosity. Nonetheless, modelinthe sinusoidal light curve, we found that, 
in both the hlters, the companion spans ~ 3.5 magnitudes between the maximum and the 
derived minimum, in agreement with previous observations. The best £t-model is shown 
as a solid curve in Figure [71 Interestingly, the light curve structure is more similar to the 
ones observed for black widow than for redback companions, which usually, b ut not always. 


show a double minimum-maximum structure due to tidal deform ation (see e.g. [Ferraro et al. 


2003al: [Cocozza et al.[ 120081: [Pallanca et al.l[2010l: [hi et al.[ [20141) . The CMD position of the 


companion during the maximum and at a mean phase (as derived by the adopted model) is 
shown in Figure [2l The system is located between the main sequence and the WD cooling 
sequence, where no normal stars are expected and thus suggesting a perturbed and strongly 
heated companion star. Again, at odds with other redback systems, this lies in a region 
more similar to t hat occupied by the two black widow companions identified so far in GCs 
(iPallancal[2014b[. Cadelano et ah, 2015). The X-ray counterpart to 47TucW shows a vari¬ 
ability which is likely due to an intra-binary shock between the PSR wind and the matter 


l ost by the companion flBogdanov et al.[[2005ll . Interestingly, as discussed bv [Bogdanov et al. 


(12006[) . the minimum of the X-ray light curve is displaced with resp ect to the optical one . 
Such a behavior has been also noticed for the black widow M71A fjCadelano et al.[ [20151 ). 
thus further strengthening the connection of this MSP with black widow systems. All this 
allow us to speculate that a scenario where 47TucW will evolve into a canonical MS P with a 
He WD companion (as in the case of MSP-A in NGC 6397; see lBurderi et al.[[2002h is some¬ 
what unlikely, opening the possibility to an evolution toward the black wido w stages. Indeed 


such a n evolutionary path has been already suggested by the simulations of [Benvenuto et al. 


(I 2 OI 4 J) . The identification of new redback companions will shed light on this possibility. 


4. SUMMARY 

By using ultra-deep, high resolution UV WFG3/HST observations of 47 Tucanae, we 
identihed the companions to four binary MSPs (47TucQ, 47TucS, 47TucT and 47TucY) and 
confirmed the two already known objects (GOM-47TucU and GOM-47TucW). The optical 
counterparts have coordinates compatible, within the errors, with the PSR nominal positions. 
In the GMD, all the objects are located in the He WD cooling sequence, as expected from the 
MSP canonical evolutionary scenario. The only exception is the companion to the redback 
system 47TucW, which is located in an anomalous region between the main sequence and 


used the “Graphical Analyzer for Time Series”, a software aimed at studying stellar variability 
phenomena, developed by Paolo Montegriffo at INAF-Osservatorio Astronomico di Bologna. 
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the WD cooling sequence, suggesting that it is a low-mass MS star highly perturbed and 
heated by the PSR flux. We compared the observed CMD positions of the detected He 
WD companions with a set of cooling tracks and derived the companion main properties (as 
masses, cooling ages, temperatures) and also some constraints on the PSR masses. All the 
companion stars have masses between ~ 0.15 Mq and ~ 0.20 Mq, and all the derived cooling 
ages are smaller than the cluster stellar population age. The orbital periods vs compa nion 
masses are in fair a greement with the evolutionary models of lTauris fc Savoniiel (119991 1 and 


Istrate et al.l (120141) . By combining the cooling age with the PSR spin down rate we found 


that the accretion history of 47TucU likely proceeded at a sub-Eddington rate. 


By taking into account our astrometric uncertainty (0.15"), we also detected a star 
having a position marginally compatible with that of 47TucH. However, its photometric 
properties would imply a PSR mass lighter than expected for a recycled NS. Moreover, its 
position in the plane of orbital period vs. companion mass is in clear disagreement with the 
theoretical predictions. While this could be due to its high eccentricity, the object could 
be just a chance coincidence and further investigations are needed before confirming its 
association to 47TucH. A possibile counterpart to 47TucI has been also identified in a star 
located in a low luminosity region of the cluster main sequence. However its distance from 
the MSP position (0.15") and the absence of any optical variability related to the orbital 
period do not allow us to asses a clear connection with the binary system. 


Finally we discussed how the properties of COM-47TucW are more similar to those 
usually observed for black widows than for redbacks, thus opening the possibility that this 
MSP could be the prototype of a redback evolving into a black widow system. 
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Table 1. 


Radio timing ephemeris of the analyzed MSPs 


MSP 

a (h m s) 

5 (° ' ") 

Offset (') 

PoRB (d) 

/ (Ms) 

Tage (Gyrs)'’ 

47TucH 

00 24 6.7014(3) 

-72 04 6.795(1) 

0.77 

2.36 

1.927x10“^ 

> 0.93 

47TucI 

00 24 7.9330(3) 

-72 04 39.669(1) 

0.29 

0.23 

1.156x10-® 

> 0.23 

47TucQ 

00 24 16.4891(4) 

-72 04 25.153(2) 

0.98 

1.19 

2.374x10-3 

> 1.43 

47TucS 

00 24 3.9779(4) 

-72 04 42.342(1) 

0.19 

1.20 

3.345x10-^ 

> 0.91 

47TucT 

00 24 8.548(2) 

-72 04 38.926(7) 

0.34 

1.13 

2.030x10-3 

> 0.32 

47TucU 

00 24 9.8351(2) 

-72 03 59.6760(9) 

0.94 

0.43 

8.532x10-^ 

2.5 

47TucW'’ 

00 24 6.059(1) 

-72 04 49.084(2) 

0.08“ 

0.13“ 

8.77x10-^ 

> 1.15 

47TucY^ 

00 24 1.4023(3) 

-72 04 41.837(1) 

0.37 

0.52“ 

1.195x10-3 

> 2.2 


Note. — From left to right: MSP name, position, offset from the GC center, orbital period, mass 
function a nd characteri s tic ag e. Numbers in parentheses are uncertainties in the last digits quoted. 
Reference: iFreire et al.l (120031 ). 


®http://www.naic.edu/~pfreire/GCpsr.html. 
’’P. Freire et al. 2015, in preparation. 
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Table 2. Optical properties of the companion stars 


Name 

a (h m s) 

5 (° ' ") 

dist (") 

niF 300 X 

mF 390 W 

COM-47TucQ 

00 24 16.489 

-72 04 25.209 

0.04 

23.19 ±0.02 

23.63 ±0.05 

COM-47TucS 

00 24 3.977 

-72 04 42.385 

0.03 

23.29 ±0.02 

23.80 ±0.05 

COM-47TucT 

00 24 8.549 

-72 04 38.965 

0.04 

23.07 ±0.02 

23.56 ±0.03 

COM-47TucU 

00 24 9.835 

-72 03 59.746 

0.06 

20.40 ±0.01 

20.85 ±0.03 

COM-47TucW 

00 24 6.063 

-72 04 49.133 

0.06 

24.28“ 

23.62“ 

COM-47TucY 

00 24 1.401 

-72 04 41.875 

0.04 

22.16 ±0.02 

22.69 ±0.04 

COM-47TucH? 

00 24 6.755 

-72 04 6.781 

0.24 

23.39 ±0.02 

24.25 ± 0.05 

COM-47TucI? 

00 24 7.953 

-72 04 39.559 

0.15 

24.14 ±0.04 

22.43 ± 0.03 


Note. — From left to right: MSP name, position, distance from the radio MSP 
nominal position, F300X and F390W magnitudes and the relatives uncertainties. 

^■The values for COM-47TucW correspond to the mean magnitudes of the best-fit 
models (see Figure [7]). 


Table 3. Derived properties of the hve MSPs with He WD companions 


Parameter 

47TucQ 

47TucS 

47TucT 

47TucU 

47TucY 

Mcom (Mq) 

~ 0.15 

n 17+0.03 
u-J.' -0.02 

0.16lHf 

0 

1 + 

0 0 
b b 

0 0 

CO to 

0.17 ±0.02 

Age (Gyrs) 

~ 5.5 

6 4 '’' 3-7 


0 .881^^ 

9 9+1.0 
^•^-1.6 

T (lO^ K) 

~ 7.6 

0 1 +1.0 

8 .0l°j 

ii.+SI 


L (10-3 Lq) 

~ 9.5 

Q 1 +1.0 

O--*-- 0.3 

10.0l°;3 

1581^7 

23.0 ±5 

MpsR (Mq) 

< 1.57 

< 4.69 

< 1.58 

< 2.30 

< 2.22 

<(°) 

> 58 

> 26 

> 57 

> 42 

> 45 



Note. — From top to bottom: companion mass, age, temperature, 
luminosity, PSR mass and inclination angle. 
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COM-47TucS 



^ 4 


« ^ 




Fig. 1.— HST images of the 2" x 2" region around the nominal position of the seven MSPs 
analyzed in this work. North is up and east is left. All the charts are obtained from a 
combinations of the available F300X images, with the exception of that of 47TucW that is 
from an image where the companion star is at its maximum luminosity. The black circles 
are centered on the radio PSR nominal position in the optical astrometric system and their 
radii are equal to our 2a astrometric uncertainty (0.30"). The red circles mark the identihed 
MSP companions. 





22 



’^F300X ^F390W 

Fig. 2.— UV CMD of the GC 47 Tucanae. Only stars with sharpness parameter |sh| < 
0.05 are plotted. The bine solid squares mark the companions to the canonical MSPs. 
The possible counterparts to 47TucH and 47TucI are plotted as an open square and circle 
respectively. Since COM-47TucW is a strongly variable object, we report its position at 
the maximum and mean luminosities, as derived by the best-£t models (see Section ITTI and 
Figure [Tj). 
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'^FSOOX ^F390W 

Fig. 3.— Same as in Figure [2l but zoomed into the WD region. The continuous curves are 
reference He WD cooling tracks for stars of 0.15 Mq, 0.17 Mq, 0.20 Mq, 0.32 Mq, 0.36 Mq 
and 0.43 Mq (from right lo left). For the two rightmost tracks, points at 1,2,4 and 9 Gyrs 
have been marked with different symbols. The photometric errors of the companion stars 
are also drawn. 
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Fig. 4.— Physical properties of 47TucQ, 47TucS, 47TucT and 47TucY (see labels), as derived 
from the comparison between the photometric characteristics of each companion and the WD 
cooling track models. In each plot, the gray lines drawn in the left and central panels show 
the allowed combinations between companion mass and cooling age or temperatnre (see 
text). The red dots correspond to the most probable valnes. In the case of COM-47TncQ 
the shaded areas mark the region (Mcom < O.ISMq) not sampled by the theoretical cooling 
tracks. In the rightmost panel of each plot, the solid curves represent the combination of 
values allowed by the PSR mass function for different inclination angles {i = 90° in red, 
i = 70° in blue, i = 50° in orange, and i = 30° in green). The blue dashed lines correspond 
to the assumed minimum NS mass (~ 1.17M0: l.lanssen et al.ll2008l ) and the largest NS mass 
value obtained for i = 90°. 
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Fig. 5.— As in Fig. 01 but for 47TucU 















26 



^COM [^©] 


Fig. 6.— MSP orbital periods plotted as a function of the best-fit companion masses, for 
each identi hed object (see lablesh plus the ones c 


etecte d in NGC 6752 and M4 (dark gray 
points; see iFerraro et ahl I2nn3al: ISignrdsson et al.l 1200311 . The three curves correspond to 
the theoretical predictions of Tauris &: Savoniie f l999h for three different stellar population 
progenitors, as reported in th e top-left leg e nd. The light gray points correspond to the 
theoretical results obtained bv llstrate et al.l (120141) . 


































0.5 1 1.5 

Orbital Phase 


Fig. 7.— Light curves of COM-47TucW in the F390W (upper panel) and F300X (lower 
panel). The two curves are folded with the radio parameters and two periods are shown for 
clarity. The black curve in each panel is the best analytical model obtained independently 
for each hlter. 




